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AN ALTERNATIVE ALGORITHM FOR COMPUTING THE BETTI 
TABLE OF A MONOMIAL IDEAL 

MARIA-LAURA TORRENTS AND MATTEO VARBARO 


Abstract. In this paper we develop a new technique to compute the Betti table of a 
monomial ideal. We present a prototype implementation of the resulting algorithm and 
we perform numerical experiments suggesting a very promising efficiency. On the way of 
describing the method, we also prove new constraints on the shape of the possible Betti 
tables of a monomial ideal. 


1. Introduction 


Since many years syzygies, and more generally free resolutions, are central in purely theo¬ 
retical aspects of algebraic geometry; more recently, after the connection between algebra and 
statistics have been initiated by Diaconis and Sturmfels in [DS98j . free resolutions have also 
become an important tool in statistics (for instance, see [Dill [SWn9| i. As a consequence, 
it is fundamental to have efficient algorithms to compute them. The usual approach uses 
Grobner bases and exploits a result of Schreyer (for more details see |Sc80l ISc91j or |Ei95[ 
Chapter 15, Section 5]). The packages for free resolutions of the most used computer algebra 
systems, like [Macaula^ Singular ICoCoA] . are based on these techniques. In this paper, 
we introduce a new algorithmic method to compute the minimal graded free resolution of 
any finitely generated graded module over a polynomial ring such that some (possibly non- 
minimal) graded free resolution is known a priori. We describe this method and we present 
the resulting algorithm in the case of monomial ideals in a polynomial ring, in which situ¬ 
ation we always have a starting nonminimal graded free resolution. A first implementation 
of the algorithm has already been tested and has led to promising results. It is our opinion 
that, by accurately refining the implementation, it is possible to achieve substantially better 
performance/timings than the already good ones shown in Tabled) 


Monomial ideals are, essentially, combinatorial objects, however the combinatorics of their 
minimal free resolutions is still mysterious. In the literature, the minimal free resolution of 
several classes of monomial ideals have been studied and fully understood. Many of these 
resolutions are cellular, i.e. they come from a (labeled) regular cell complex built by the 
generators of the monomial ideal (see |BS98j or [MSOdl Chapter 4]). Due to this reason, 
cellular resolutions have been intensively studied in the recent literature, giving rise to a 
new interesting area of research. However, not every monomial ideal admits a minimal 
cellular resolution, though all monomial ideals admit some (not necessarily minimal) cellular 
resolution. In this paper, we use algebraic discrete Morse theory to algorithmically reduce 
a cellular resolution of a given monomial ideal to its minimal free resolution. We explicitly 
describe the algorithm when the starting cellular resolution is the free resolution found by 

1 



























2 


M. TORRENTE AND M. VARBARO 


Lyubeznik in |Ly88| , which is actually simplicial (i.e. the regular cell complex is simplicial). 
Our implementation of the algorithm is provided in this case. 

On the way of writing the algorithm, we prove some rigidity properties, like the following 
vanishing for monomial ideals generated in degrees < d fTheorem 16.21) : 

Pi,k{S/I) = 0 V fe = j,..., J + d - 1 ^ Pi+ij+d{S/I) = 0. 

This strengthens the recent result proved by Herzog and Srinivasan in [HS14j saying that, if 
l3i,k{S/I) = 0 V /c > j, then Pi^i^k+diS/1) = 0 \/ k > j. Theorem 16.21 has independently been 
proved, by different techniques, in |Yal5j . 

Finally, a similar idea could be used to compute the simplicial homology with coefficients 
in a field of a finite simplicial complex (see Remark 14.51 for more details). Already many 
algorithms are available for this task; nevertheless, we think it would be worth to inquire on 
the potentiality of this idea. 


2. Reducing resolutions 


Let i? be a ring, and F, = 
natural numbers {0,1, 2,...}. 


a complex of free R-modules, where N is the set of 
For any i G N let A* be a set to index the basis of F), that is: 


Fi = 0 R • CT. 

(tEA* 

With respect to such a basis, we can write the differentials dj : Fj — )■ -Pj-i as: 

Icr ^ [o- : r] • Ir, 

rGAj_i 

where [a : r] is an element of the ring R. 

We are going to describe a special case of the so-called Algebraic Discrete Morse Theory, 
an useful tool to reduce the size of F, (when possible). Let a G Aj and (d G Aj_i such that 
[a : /?] is an invertible element in R. Then construct the following complex F, = (Fj,5j)jgN 
of free R-modules: 


Fi = 


© 

© 

Fi 


aGA, \{«} ^ ^ 
reA,_i\{/3} R ■ ^ 


if i = j, 
if i=j -1, 
otherwise 


The differentials di : Fi ^ R-i are defined as: 


di{la) 


di{la) 

< Z:rGAj\{o}[^ : t] • Ir 
ErGA,_i\{/3} (k : 


[a-.p] ) 


if ii{j + l,j}, 
ifi = j + l, 
ifi = j 


The following lemma is a particular case of [.TWO91 Theorem 2.2]. 


Lemma 2.1. F, = (Rj,5j)jgN is a complex of free R-modules homotopically equivalent to 
F, = (Ri,9j)jgp}. In particular, 


Ri(F.) = Ri(F.) ViGN. 
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Let US remind that F, is a free resolution of an i2-module M if Hi(F,) = 0 for any i > 0 
and Hq{Y») = M. Therefore the above lemma implies that, if F, is a free resolution of M, 
then F, is a free resolution of M as well. 

We are interested in the case that R is a polynomial ring over a field k and F, is a graded 
free resolution of a finitely generated graded i2-module M. In this case, F, is not minimal 
if and only if there exist a G Aj and (5 G Aj_i such that [a : /3] G k^. Therefore Lemma 
12.11 provides an iterative procedure to get the minimal graded free resolution of M if some 
graded free resolution of M is known a priori. 


3. SiMPLICIAL RESOLUTIONS 


Let S = k[a:i,..., x„] be the polynomial ring in n variables over a field k, and / C ^ be 
a monomial ideal generated by monomials u = ui,... ,Ur. Let A be a (labeled) simplicial 
complex on the set {1,... , r} =: [r], where the label of a face a G A is the monomial rria G S, 
where: 


(1) rria := lcm(ttj|i G <t) G S'. 

Let us remind that the reduced simplicial homology of A (with coefficients in k) is the 
homology of the complex of k-vector spaces C,(A; k) = (Cj, where d — 1 is the 

dimension of A, 

Ci = k • a 

a 

where the sum runs over the ^dimensional faces of A, and the differentials are defined as: 

1(7 ^sign(u,cj) • 

vdcr 

where 


( 2 ) 


sign(u,cr) := ( — 1)'^ ^ G k 


if V is the qth element of cr. By imitating this construction, let us consider the graded complex 
of free 5-modules F,(A;;u) = (F), where 

Fi = 05(-deg(m(^)) 


where the sum runs over the {i — l)-dimensional faces of A, and the differentials are defined 
as: 


1(7 


E 


. , . rn„ 

sign(u, fj j- 

^a\{v} 


• l(7\{w}- 


Since any v G [r] is a vertex of A, one has that Hq(F,{A]u)) = S/I. Furthermore, because 
the differential are graded of degree 0, one can check that F,(A;m) is a graded free resolution 
of S/I if and only if A<u is an acyclic simplicial complex for any monomial u G 5, where 
A<u consists in the faces <7 G A such that ma\u (cf. |BPS981 Lemma 2.2]). We will call this 
kind of resolutions, introduced for the first time by Bayer, Peeva and Sturmfels in |BPS98| . 
simplicial. In particular, if A = /C is the (r — l)-simplex, then F,(/C;m) is a graded free 
resolution of S/I: indeed in this case A<„ is the simplex on the vertices i such that Uj|u, 
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which is contractible. That F,(/C; u) was a resolution was proved, for the first time, by Taylor 
in |Ta60] . 

Definition 3.1. With the above notation, F,(/C;;u) is called the Taylor resolution. 

The Taylor resolution is far to be minimal in general: in fact, it tends to be huge. A 
smaller graded free resolution, still simplicial, was found by Lyubeznik in |Ly88| . For each 
k G [r] and u C [r], denote by 

(3) (T>fc := {i ^ a : i > k}. 

Let C be the simplicial complex on [r] consisting of those faces a C [r] such that does not 
divide for any k & [r]. Then Lyubeznik proved in |Ly88| that F,(£;u) is a graded free 

resolution of S/I. 

Definition 3.2. With the above notation, F,(£;;u) is called the Lyubeznik resolution. 

Notice that the Lyubeznik resolution depends on the way in which we order the monomials 
ui,... ,Ur. Even if the Lyubeznik resolution tends to be much smaller than the Taylor 
resolution, it may happen (indeed it often happens) that it is not minimal for any order of 
the generators of I. 


4. The algorithm 

Let I T S = k[xi,..., Xr/\ be a monomial ideal generated by monomials of positive degree 
wi,..., Ur- (r > 0). If cr = {ii,..., is} and r = {ji,..., jg} are subsets of [r] of the same 
cardinality, by cr < r we mean that a is less than r lexicographically, i.e. that there is k such 
that ik < jk and ip = jp for all p < k. For any i G [r], define the following set: 

(4) Lyuj := {cr C [r] : \a\ = i and 3 A: G [r] : Uk\m^^^}, 

where cr>fc was defined in ([3]). Construct the graph G{I) as follows: 

(i) The set of vertices of G{I) is Aq U Ai where 

Aj := (cr C [r] : |cr| = z, cr ^ LyUj} V z = 0,1. 

(ii) There are no edges. 

The purpose of the following algorithm is to iteratively transform the weighted graph G{I) in 
a way that will provide the Betti table of S/I (see sections [2] and [3]) . During the algorithm, 
the vertices of G{I) will be always partitioned in 2 sets, namely Aj and Aj+i. The main 
loop of the algorithm works as follows. If either Aj or Aj+i is empty, the algorithm stops. 
Otherwise, the set of vertices Aj is stored and deleted, and a new set of vertices Aj +2 is 
constructed together with a set of weighted edges A*^^, going from Aj +2 to Aj+i. Both sets 
Aj _|_2 and A*^^ continually change, and eventually the edges of A*^^ will have weight 0. At 
the end of each loop, the counter i is increased by one. In the worst case, when z = min{n, r}, 
the set Aj+i is empty, so the algorithm stops. 
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Remark 4.1. With the above notation, the following inequalities follow by the Taylor’s 
resolution: 

projdim(S'//) < min{r, n} and reg(£'//) < projdim(£'//) • max{deg(Ri) — 1} 

i 

In particular, the Betti table of S/I has size at most a{I) x h{I), where a{I) := min{r, n} + 1 
and b{I) := min{r, n} • maxj{deg(uj) — 1} + 1. 

Algorithm 4.2. Let Aq, Ai, a(/), 6(7) be defined as above; set i := 0, Vi := Ai, and let B 
be the (a(/) x 6(/))-matrix with all zero entries. 

SI: Let Bi be the (o(I) x 6(/))-matrix whose entries are all zero except for the (i + l)-th 
column: the (j, i + l)-th entry will be 

\{a £ Ai : deg(m^) =i + j- 1}|. 

Analogously, let Rj+i be the (a(/) x 6(/))-matrix whose entries are all zero except for 
the {i + 2)-th column: the (j, i + 2)-th entry will be 

|{cj G Aj+i : deg(mcr) = i + j}\- 

S2: If i = r — 1 then set B := B + Rj + Rj+i and stop; 
else if i = n then set B := B + Rj and stop; 
else if Aj = 0 then stop; 

else if = 0 then set B := B + Rj and stop; 

else set B := B + Rj and remove the set of vertices Aj from G{I). Change the 
weighted graph G{I) as described in the procedure The deformation of G{I) below, 
set i := i + 1 and go to step SI. 

The deformation of G{I) 

Set 14+2 := {t U {u} : r G 14 + 1 , £ H \ t}, Aj +2 := 0 and A*^^ := 0. Iteratively transform 

Aj +2 and A*^^ by lexicographically running on the elements a G 14+2 in increasing order as 
follows: 

• If cj G LyUj+ 2 , then I 4+2 := 14+2 \ {<t} and go to the next a. 

• Otherwise, construct the following sets: 

(Af+i)i := {a \ {u} G A^+i : u G a, and (u, a \ {u}) ^ A'+^j 

(A-+i) 2 := {/3GA,+i:(a,/3)GA*+?} 

(A-+i) 3 := {/3 G (Af+i)2 : [a :/3] ^ 0} 

• If := (A^_^;^)i U (Af _^_;^)3 = 0, then set Aj +2 := Ai +2 U {a} and go to the next a. 

• Otherwise, for every j3 G (A^_,_;^)i, set [a : /3] := sign(iT \ (3,(3) (see ([2])). 

Let r be the smallest element of A^^^. Reset A^^_^ := A^^^ \ {r}. 

Set A^+^ := (A^+^)i U (A(+^) 3 , where: 

(A(.+^)i := {a = r U {u} : u G [r] \ T, a > a, Uy\m„, a ^ LyUj+ 2 } 

(A;+2)2 := {a G Ai +2 : a > u, (a, r) G A^+^j 
(A;+2)3 := {aG(A;+2)2:[«,r]^0} 
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If a G (A^+^)i \ (A^+^) 2 , define [a : r] := sign(a \ r, a). 

Set Aj +2 := Aj +2 U (A^+^)i and, for all (a,/3) G A^+^ x precede as follows: 

- If (a,/3) G A*+2^ then [a,/3] := [a : /3] - ■ 

- If (a,/3) i A*+2, then A*+2 := A*+2 y {(a,/3)} and: 

* If a D /3, then [a, /3] := sign(a \ /3, a) - 

* Otherwise, [a,/3] := — ■ 

Set Ai +2 := Ai +2 \ ({o'} n Ai+ 2 ) and Aj+i := Aj+i \ {r}. Put 

A*+? := A*+2 \ ({(«,^) : « e (A;+2)2} u {(a,/3) : /3 G {^^M) 
and go to the next a. 


Remark 4.3. Note that the above algorithm can be modified to provide the multigraded 
Betti numbers of 5/1: namely, given a vector a = (ai,..., a^) G Z"', 

/3i,a{S/I) = \{a e Ai : rua = x1^ ■ ■ ■ xl"}\ V i G N, 

where Aj is taken at the step i. 

From the algorithm described above one can immediately infer some interesting upper/lower 
bounds for the depth/Castelnuovo-Mumford regularity of S/I. With the notation above, let 
I = (til, ■ ■ ■) Ur) C S' = lk[xi,..., Xn]. We say that a subset a C [r] is critical (for I ) if: 

mo- / rnrr\{v} ^ v ^ a and mo- ^ mou{^} V u G [r] \ cr. 

We introduce the following two invariants of a monomial ideal I: 

(5) p{S/I) := max{|(T| : a is critical} 

(6) r{S/I) := max{deg(mo-) — \a\:a is critical} 

Proposition 4.4. If a is a critical set for I, then l3\(T\,deg(ma){S/1) 0. In particular: 

depth(S//) < n — p{S/I) and reg(S//) > r{S/I). 

As a special case, we recover a result of Katzman [Kan6| : let G be a simple graph of 
n vertices and I = {xiXj : {t,j} is an edge of G). He noticed that, if G has an induced 
subgraph H consisting of t disjoint edges, then reg(S//) > t. Clearly, in such a situation the 
set a corresponding to the edges in H is critical for I, so that r{S/I) > deg(mo-) — |(t| = t. 

Remark 4.5. To compute the multigraded Betti numbers of a monomial ideal, we can always 
reduce ourselves to square-free monomial ideals by polarization. A square-free monomial ideal 
/ = {ui, ... ,Ur) C S is associated to a simplicial complex A on n vertices via the Stanley- 
Reisner correspondence. Thanks to the Hochster’s formula (cf. [MSOdl Corollary 5.12]), to 
know the multigraded Betti numbers of S/I is equivalent to know the simplicial homology 
with coefficients in k of all the 2” induced subcomplexes of A. 

If one is only interested in the simplicial homology of A, it is enough to know the multi¬ 
graded Betti numbers A,(i,...,i)(5'/-I), more precisely: 

dimkRi(A;k) = Pn-i-i,{i,...,i)iS/I) 
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Being our algorithm “local”, to know such multigraded Betti numbers is not necessary to 
compute the minimal free resolution of S/I completely. In fact, it is possible to implement 
a faster version of the algorithm aimed at the computation of the simplicial homology with 
coefficients in a field. Already several algorithms are available for this task, with applications, 
for example, in network theory. Nevertheless, we wish to study the potentiality of such an 
algorithm in a forthcoming paper. 

To give a preview of how such an algorithm would work, let {Gi, ..., G^} be the set of 
minimal nonfaces of A. Essentially, one should perform the “deformation of G(/)” only on 
the sets of vertices 

Vi = {c7 C [r] \ LyUj :\a\ = i and [_J Gj = [n]} i = n — dim A,..., n. 

jGcr 

Then dimt fIi(A; k) will be the cardinality of the set A„_i_i (with the same notation of 
the algorithm) at the end of the deformation. A remarkable feature is that to compute the 
simplicial homology in this way, one needs as input the minimal nonfaces of the simplicial 
complex, rather than its facets. This fact is convenient in certain situations. For example, 
in many cases is of interest to compute the homology of the clique complex of a given graph: 
while the facets of such a simplicial complex are the maximal cliques of the graph, its minimal 
nonfaces are simply the pair of nonadjacent vertices, which are easier to store if only the graph 
is given (as it typically happens in network theory). 

5. Experiments 

In this section we present some examples to show the effectiveness of our algorithm. We 
implemented it using the C++ language and some routines available in CoCoALib, a GPL 
C++ library which is the mathematical kernel for the computer algebra system CoCoA-5 
|CoCoAj . All computations are performed on an Intel Core 2 i5 processor (at 1.4 GHz) 
using both our new algorithm and the function BettiDiagram included in CoCoA-4. We 
are not comparing our timings with CoCoA-5 itself because the implementation for minimal 
free resolutions has not yet been optimized in the current version CoCoA-5.1.1. We perform 
three experimental tests: in the first case we consider monomial ideals of degree d = 2 (results 
summarized in Table l.(a)); in the second case we consider monomial ideals of degree d = 4 
(results summarized in Table l.(b)); in the third case we consider monomial ideals of degree 
d, with 5 < d < 8 (results summarized in Table l.(c)). 

Tables l.(a), l.(b) and l.(c) consist of four columns: the first column, labeled with n, 
contains the number of variables of the polynomial ring; the second column, labeled with r, 
contains the number of minimal generators of the monomial ideal; the last two columns, 
labeled with “New” and “CoCoA”, contain the average time to compute the Betti tables 
of 10 monomial ideals using our new algorithm and the software CoCoA-4 respectively. In 
all tested numerical experiments our algorithm reveals to be more efficient than CoCoA- 
4, since it returns the Betti tables in a shorter computational time. We also performed 
some comparisons of our results with computations done using the software Macaulay-2 
(see |Macaulay2| ), in which case the computational timings are comparable, with variations 
depending on the size of the example. In our opinion, though, by accurately refining the 
implementation, the efficiency of our algorithm has substantial room for improvement. 






(a) Monomial ideals of degree 
d = 2. 



Timings 

n 

r 

New 

CoCoA 

10 

8 

0.001 

0.02 

10 

10 

0.005 

0.04 

10 

12 

0.014 

0.11 

15 

12 

0.023 

0.19 

15 

15 

0.177 

1.32 

15 

18 

0.512 

2.33 

18 

18 

2.111 

13.14 

18 

20 

9.608 

32.87 
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(b) Monomial ideals of degree 
d = 4. 



Timings 

n 

r 

New 

CoCoA 

10 

8 

0.002 

0.02 

10 

10 

0.004 

0.05 

10 

12 

0.012 

0.11 

15 

12 

0.029 

0.43 

15 

15 

0.208 

2.91 

15 

18 

0.786 

9.37 

18 

18 

3.282 

78.88 

18 

20 

24.194 

158.27 


(c) Monomial ideals of degree 
d, with 5 < d < 8. 




Timings 

n 

r 

New 

CoCoA 

10 

8 

0.002 

0.04 

10 

10 

0.007 

0.09 

10 

12 

0.014 

0.23 

15 

12 

0.035 

0.81 

15 

15 

0.261 

7.51 

15 

18 

1.247 

58.84 

18 

18 

4.064 

271.5 

18 

20 

63.00 

482. 73^ 


Table 1. Tests on monomial ideals of degree d = 2, d = 4 and 5 < d < 8 


6. A VANISHING INSPIRED BY THE ALGORITHM 


In the last section of the paper we prove a vanishing result arisen while developing the 
algorithm. Before stating it, for the convenience of the reader we recall how to get the 
Mayer- Victories spectral sequence for simplicial homology: 


Remark 6.1. Let A be a simplicial complex on [n], and Aq, ..., A^ be subcomplexes such 
that 

A = U A.. 

i=0 

For any sequence of integers 0 < qq < • • • < Rp < •s, let us denote by 


p 



k=0 


Fixed a ring R, for any simplicial complex F we denote by C,(F; R) the complex of R-modules 
defined as in Section [3] in the particular case that R was the field k, and by C*(F;R) := 
Homj^(C,(F; i?)). One can check that the following is an exact sequence of complexes of 
ii-modules: 


0 ^ C*(A; R) ^ C‘’°{R) ^ C‘’\R) 


d2 


^C'’^-\R) ^0, 


where C‘’P{R) = ®ao<...<ap C*(Aao,...,ap; -R) and, for a given j, if a = (Q;ao....,aj,)ao<...<aj, is an 
element of C^'P{R), the {bo,... ,6p+i)th component of dP~^^{a) € C^’P^^{R) is: 


P+i 




A:=0 


a. 


)r v?^p+l 


C'-’ (^6o,...,6p+i ) 


^this average is made on 7 tests out of 10, since in 3 experiments the computation ran out of memory 
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(considering A as A 0 ). Let us consider the double complex of i?-modules where the 

differentials 

for a given natural number q, are the ones used to define the simplicial cohomology of ^ao,...,aq 
with coefficients in R. By the general theory of spectral sequences arising from double com¬ 
plexes (c.f. [GM031 Chapter III.7, Proposition 10]), we get the following spectral sequence: 

(7) 0 ^ 

0<ao<...<ap<s 

Given two elements a = (oi,..., a„) and h = ( 61 ,..., 6 „) of Z"", for a < b we mean a* < 5j 
for all i €z [n]. Furthermore, for a < b we mean a <b and a ^ b. 

Theorem 6.2. Let I = (ui,..., Ur) be a monomial ideal of S = k[a;i,..., Xn], and let = 
degin{u() for all £ £ [r]. For any ^ € [r], i £ N and a £ Z"; 

Pi,b{S/I) = 0 ^a<b<a + ai /3i+i,a+a^(5'//) = 0 . 

In particular, if d = max£{degg(u^)}, for any j £ Z.- 

(3iMS/I)=0 ^ k = j,...,j + d-l A+i,i+d(5//) = 0. 

Proof. It is not difficult to see that, by polarizing, one can reduce himself to the square-free 
case (though one has to be a bit careful because by polarizing the multi-degree changes). In 
such a situation, for any s £ N, j3s,c{S/I) / 0 only if c = (ci,..., Cn) is a square-free vector, 
that is Cj £ {0,1} for any j £ [nj. It is convenient to think at such square-free vectors as 
subsets of [n], namely 

c ^cr(c) = {j £ [n] : Cj = 1 }. 

Notice that c\ < C 2 if and only if cr(ci) C cr(c 2 ). Let a, ae and 7 = a U be the subsets 
of [n] corresponding to, respectively, a, and (a -|- a£)(mod2). Furthermore, let A be the 
simplicial complex such that I = I/s. and set F = A..^. Notice that our vanishing hypothesis 
can be written as: 

l3i,^{S/I) = 0 y aC 13 Cj. 

By Hochster’s formula, this is equivalent to say: 

( 8 ) iLl^l"*“^(r^;k) = 0 V a C/3 C 7 . 

Since is not a face of A, 

U ~ 7". 

Obviously, for all subsets 6 of af. 

vGS 

Therefore the Mayer-Vietoris spectral sequence d?]) with respect to the covering (F..|,\{^})ugc^ 
of F is: 

SGa£ 

\s\=p+l 
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So, by ([8]) we infer that * ^(r;k) = 0. We conclude by using Hochster’s formula that 


A+1,^(5//) = 0. 


□ 


Remark 6.3. As a particular case of the above theorem, we get: 

= 0 V fc > j ^ A+i,fc+d(5//) = 0 V fc > j. 

This fact has recently been proved by Herzog and Srinivasan in [HS141 Corollary 4]. A further 

consequence is the following result of Fernandez-Ramos and Gimenez [FG141 Theorem 2.1]: 

if I is the edge ideal of a graph, then j3ij(S/I) = /3i ,,+1(5/1) = 0 ^ A+i,,+2(5/1) = 0. 
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